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Abstract 



This paper is an investigation of the effects of a thermal bulk fluid in 
brane world models compactified on AdS^. Our primary purpose is to study 
how such a fluid changes the bulk dynamics and to compare these effects 
with those generated by matter localized to the branes. We find an exact 
cosmological solution for a thermally excited massless bulk field, as well as 
perturbative solutions with matter on the brane and in the bulk. We then 
perturb around these solutions to find solutions for a massive bulk mode in 
the limit where the bulk mass (m^) is small compared to the AdS curvature 
scale and T < rriB- We find that without a stabilizing potential there are 
no physical solutions for a thermal bulk fluid. We then include a stabilizing 
potential and calculate the shift in the radion as well as the time dependence 
of the weak scale as a function of the bulk mass. It is shown that, as opposed 
to a brane fluid, the bulk fluid contribution to the bulk dynamics is controlled 
by the bulk mass. 
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I. INTRODUCTION AND MOTIVATION 



The possibility that we are hving on an isolated defect has led to new paradigms which 
are now being explored. While the low energy effective theory of these models correctly 
reproduces the standard model physics, we expect that the cosmology should be quite dis- 
similar at large enough temperatures. In particular, the cosmology of models which solve the 
hierarchy problem through warped compactification, such as the Randall Sundrum model 

(RSI), are especially interesting because we expect to see deviations from the standard 
cosmology as we approach the weak scale [0-0. The cosmology of these models has been 
discussed previously including the effects of matter density on the branes. What makes the 
cosmology of these models particularly interesting is that we expect five dimensional gravity 
to become strongly coupled at the TeV scale. Above the TeV scale, it is believed that the 
theory should still have a description in terms of a local quantum field theory, via the holo- 
graphic principle In this paper we will continue the exploration of the cosmology of 
(RSI) by studying the effects of a thermal fluid in the bulk. We will study temperatures low 
enough so that, T < nikkM^/k (where M5 is the five dimensional Planck scale and k is the 
AdS curvature scale) and a weakly coupled {k < M5) five dimensional gravity description 
is still valid. Since our results only depend upon the geometry they are also applicable to 
models where the hierarchy is stabilized via supersymmetry and generated via the geometry 
i TSj . However, we will couch our results in terms of (RSI). 

The primary motivation for this work stems from the fact that there must be some 
nontrivial bulk physics in the RSI model, as well as in the supersymmetric models mentioned 
above. Thus, we would expect that at some point in its thermal history, more specifically at 
the TeV scale, a thermal bulk will be relevant. In its purest form, RSI entails the trapping 
of the standard model particles on a brane of negative tension positioned at an orbifold 
fixed point. At some finite proper distance in the extra dimension there is a positive tension 
brane where the graviton is localized. In this choice of coordinates, the small overlap of 
the graviton wave function with the standard model brane leads to the apparent weakness 
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of gravity. This minimal model must be augmented to account for the stabilization of 
configuration. That is, in the minimal RSI, the distance between the branes parameterizes 
a flat direction in field space. The introduction of new physics in the bulk can lift this fiat 
direction. For instance, the Goldberger-Wise (GW) mechanism is a natural way to lift the 
degeneracy by introducing a bulk scalar with a non-trivial vacuum profile. Thus, in general, 
we may expect that in addition to gravitons propagating in the bulk, there will be other 
fields which may play a role in the bulk physics. In fact, it is phenomenologically viable to 
have particles with standard model gauge charges propagate in the extra dimension ||16|. If 



this were the case then bulk thermal effects would be enhanced due to the large number of 
species which would equilibrate. 

In this work we seek to understand how the cosmology changes once the bulk thermal 
effects are taken into account. We would expect that at long distances the bulk fluid will 
have the same effect as a fluid trapped to the brane. Though we might be concerned by 
the fact that, since the lapse function is changing exponentially, the effect of the bulk fluid 
at long distances could differ from those of a brane fluid. It has been pointed out that the 
visible and Planck brane matter densities (p, p*) contribute to the expansion rate as ^ 

H- = ^-^{p. + pe-^''% (1) 

where is the proper distance between the branes. Given that the natural scale of the 
Planck brane is M^;, this puts strong constraints on the brane matter density, as pointed 
out in [0 . Indeed, this constraint implies that the Planck brane must be in its vacuum state 
for all intents and purposes. Thus, depending on the bulk matter density profile, a bulk 
energy density could lead to over-closing the universe. 

We would also expect a bulk fluid to effect the bulk dynamics quite differently from a 
brane fluid. A brane fluid generates a non-trivial radion potential, which in the absence of 
a stabilizing mechanism, tends to push the branes apart, unless one imposes the unphysical 
fine-tuning between the visible brane energy density (p) and Planck brane density (p*), 
p^ = — pe~^'^''° [0^- Here we calculate the radion potential generated by a bulk thermal 



fluid |T8[, and determine its effect on the above fine tuning relation. We then include a 
stabilizing potential to see how a bulk fluid shifts the radion expectation value and the weak 
scale as functions of time. 

To address these issues we will solve the coupled bulk fluid Einstein equations. We start 
by finding an exact solution with fluid in the bulk and vacuum branes. In general, this is 
rather intractable problem, since the sources on the right hand side of the Einstein equation 
are complicated highly non-linear functions of the metric ansatz. However, we will be able 
to find an exact solution in the particularly simple case of a massless field at temperatures 
too low to excite higher Kaluza-Klein modes. We then find a realistic solution with matter 
in the bulk, as well as on the branes, which is valid when the energy densities are small 
compared to the bulk curvature (i.e. visible brane temperatures less than a TeV). We then 
perturb around the massless solution to find more interesting solutions which include the 
effects of a massive bulk field. We show that without a stabilizing mechanism, there are no 
physical solutions. Furthermore, once stabilization is taken into account, all the non-trivial 
bulk dynamics, e.g. the evolution of the weak scale, are driven by the bulk field mass. 



II. EXACT SOLUTION WITH A BULK FLUID 

The five dimensional action is given by 

L = J d'xy^{-MlR -A + Lb) + J <fx^ LreV + Lpi, (2) 

where g and are the induced metrics on the visible and Planck branes, respectively. M5 
is the five dimensional Planck scale and Lb is the bulk field Lagrangian which describes 
the stress energy of the bulk fluid. The fifth dimension is compactified on S1/Z2, with the 
Planck and visible branes placed at the orbifold fixed points, y = 0,1, respectively. For this 
section we will leave the branes empty, and thus, their only contribution to the action are 
due to their tensions and V. 

We will search for solutions which have a stable (time independent) fifth dimension. We 
thus make the following ansatz for the metric. 
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We could rescale the y coordinate to eliminate the constant 60, but choose not to for book- 
keeping purposes. The bulk stress energy tensor is 

With TP = diag(pB(t,y), -pB{t,y), -P^{t,y)). 
The Bianchi identity leads to 

A'^{y){pB{t,y) -ZPB{t,y) + AP^{t,y)) - Pl,{t,y) = 0. (5) 

PB{t, y) + 3-(t) {PB{t, y) + pB{t, y)) = 0. (6) 
a 

Primes denote derivative with respect to and dots denote derivative with respect to the 
time-like coordinate. 

The Goo, Gii and (^55 Einstein equations in the bulk are 

(^)'-^e-^°^^)(2K(y)f =^exp-^°(^)(p^^ , (7) 

^^e-^My)i2{A^^y)f_A^^y))^{0^^ , (8) 

6(^o(yyr - 36^e^^o(^) {0{t) + l{t)^ ^ KX{P.{t,y) - K) , (9) 
and we have defined n"^ — Combining the Goo and Ga equations, we see that the bulk 
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energy density and pressure must be of the highly restricted form 

PB{t, y) = p{ty^'^'\ P^it, y) = P,{t)e^^'^^^ PB{t, y) = PB{t)e^^'^^\ (10) 

This restriction results from our simple choice of ansatz. Whether or not a bulk fiuid will 
indeed yield this form will be discussed later. 

To solve this system of equations wc first make linear ansatz for Aq, since the geometry 
should reduce to AdS in the limit of vanishing bulk energy-momentum, which leads to the 
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usual result []Ao = ^&o|l/| with k = — as in the RSI solution. This result has been 
derived using the Goo and Ga equations. We have not, to this point, imposed the G55 
equation which for this simple ansatz is identical to the Bianchi identity 

2p,it) + psit) - 3pB{t) = 0, (11) 

With the result that 



pBit) = -A^m^ (12) 



PB{t) = —A[-m+2-{t), (13) 



k'^ Way ' ' a 



which is seen to obey (^. A related solution was found previously in pil|JTT|] . Note that the 
fluid tends to build up near the visible brane. It is straightforward to see that the usual 
cosmology on the brane results from this solution. 

We must now deal with the fact that our ansatz leads to the constraint ([lOD on the form 
of the bulk pressures and energy density. One simple case for which we will have a solution 
is vacuum domination, as discussed in [|10| , pT| , where the constraints (|T0|) are obviously 



satisfied. Here we are more interested in the more inevitable case where the energy density 
arises from thermal fluctuations. In this case we must calculate the local pressure and energy 
density using our solution to determine if the constraint is indeed satisfied. We will postpone 
this calculation until we have included the effects of matter on the branes. 
Before closing this section we note that the jump conditions 

3[A'M]- = ^'V. - 3[A;(1)]- = -«:V, (15) 



From here on, the absolute value, implied by the orbifolding, will be dropped for notational 
simplicity. 
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lead to the usual fine tunings. One tuning is needed to fix the flatness of the brane and the 
other to enforce the zero force condition between the branes, 

K = Qk/K\ V = -K. (16) 

The fact that there are two fine-tunings is a consequence of the lack of a radion potential 
for this solution, which one might have expected to be generated by the existence of the 
bulk fluid. This result is clearly related to the constraint ([10|) . We shall explore this point 
further below. 



III. SOLUTIONS WITH MATTER ON THE BRANE AND IN THE BULK 

Solutions with matter on the brane were previously found in Here we augment 

these results to include the effects of a massless bulk fluid. The derivation is similar to that 
in but is included for completeness. 

The previous ansatz is not general enough to allow for a solution with brane matter. We 
thus introduce the brane stress energy tensors 

Tb = ^diag(K + P., K - P., K - P.. V. - P., 0) + ^^^^diag(- V + p, - \/ - p, V - 0), 

(17) 

and we modify the ansatz, as follows 

ds^ = exp-2^»(^) (1 + 2 fit, y)) - a{tf exp'^^^^^) (1 + 2g{t, y)) dx" - bldy^ (18) 

and treat both the bulk and brane matter as perturbing sources around the standard RSI 
solution. We will calculate to leading non-trivial order in the perturbations p ~ 0{e). 
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In [12 1 a solution was found which is accurate to all orders in the brane matter density. However, 
the solution is only valid as long as the energy density on the brane is small enough that it can 
still be assumed that the G55 Einstein equation is automatically satisfied. This will be true only as 
long as the matter energy density is small compared to the scales in the stabilizing field potential. 



Since p ~ ^p, it will often be the case that we will be able to drop time derivatives of the 
fluctuations. The solutions will again restrict the bulk profile of the matter density to be of 
the form ([10|). 

The Go5 Einstein equation is now no longer automatically satisfied, and is given by 

f\t,y) + -{t){f'{t,y)-g'{t,y)) = 0. (19) 
ct 

In this equation we must keep the time derivative of the perturbation, since it is leading 
order. While the (00), {ii) and (55) linearized Einstein equations are given by 

3 {l)'^'^"^'^ - HAM? + l2AMf'it,y) - 3f"it,y) = by{e-'^^^y^pB{t,y)), (20) 
- 4^(2/'(t, y) + g'{t, y)) + ^{V"{t, v) + /(t, v)) ' ^^it) + (^)(t)) = n'{Ps{t, y)), 
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SbyMy) (f^\t) + -{t)]-3AMi3f\t,y)+g'{t,y)) = blK'P'^{y,t). (22) 



The y component of the contracted Bianchi identities is still of the form (^ and the energy 
conservation equations for the brane as well as bulk matter take on the form of (^. 
The jump conditions are 



[/'(O)]- = -^P.(t), [/'(I)]- = ^p(t). (23) 



[g'm. = + \p.it)). [9\^)]- = -^%P{t) + \p{t)). (24) 

We may now solve for f{t,y) and g{t,y) as follows. As before the leading order Einstein 
equations lead to Ao^y) = kh^y. Following our exact solution in the previous section, we 
assume the matter has the form (plQl). Then integrating the next to leading order Goo 
equation leads to 
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f{t^y) = A{t)e^''y 



a) 6 



e^^'y. (25) 



Imposing the jump conditions at and 1, yields the relations 



8k^A{t) 



(26) 



Pit), 



(27) 



which leads to the FRW-hke equation for the expansion rate 

2 StiG 



(28) 



Where we have used the relation SttGn — i_'^^2kbo , which can be read off of the effective four 
dimensional action. To leading order, ^ is the expansion rate on the visible brane. Notice 
that this result agrees with the contribution to the expansion rate that one would get by 
averaging the bulk energy over the extra dimension. As with the Planck brane matter there 
is no warp factor suppression, which may at first be alarming. To determine whether or not 
this result presents a problem in these models we need to know how the bulk energy density 
is normahzed. Obviously, if it scales hke Mpi, the cosmology is not viable. The correct 
normalization will be calculated in section four. 

We may now solve for g{t,y) by assuming an equation of state ps — wbPb and imposing 
the Go5 equation in conjunction with the time component of the Bianchi identity. 



-9{t,y) 



^2kboy 



4:k^ 



d 
~dt 



+ e^'=^°ni(t) + -A(t)) 
a 

(29) 



Imposing the jump conditions for g{y)i reproduces the usual energy conservations for the 
branes 



p{t)+^-{t){m+p{t)) =0 



pS) + ^-it) {P.{t) + p.{t)) = o. 

a 



(30) 



Thus, to the order we are working, the branes expand at the same rate. 

Similarly the Ga Einstein equation leads to the FRW like relation for the acceleration 
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K K 

(31) 

The G55 equation leads to the requirement 

AA{t) + -A{t) = 0. (32) 
a 

However, upon inspection we see that A{t) is independent of the bulk energy density and 
pressure and that the constraint from the G55 equation is just the usual unphysical relation, 
which results from the lack of a stabilizing potential, p*(t) = —e~'^^^°p{t). This is just 
a sign that the bulk fluid has generated no radion potential. We shall see that this is a 
consequence of our simple ansatz for the metric, which greatly restricts the form of the 
bulk energy density, and that this solution will only be valid for a massless bulk field at 
temperatures low enough as not to excite higher KK modes. 

IV. THE BULK ENERGY DENSITY 

The solution found in the previous section necessitated that the bulk energy density and 
pressure have the particular form (0). As we mentioned previously, while the constraint is 
satisfied for a bulk cosmological constant we are more interested in the case of a thermal 
fiuid, as this is the case of phenomenological relevance. We thus would like to calculate the 
thermal expectation value 

{TABit, y)) = Tr[TAB{t, y)e-*^], (33) 

in our background solution to test for consistency. We have assumed that the expansion 
rate is small enough that we have an approximate time-like killing vector, and that local 
thermodynamic equilibrium can be established. Furthermore, we normalize the time-like 
killing vector such that the Hamiltonian generates time translations on the visible brane. 
That is to say that, T will correspond to the proper temperature on the brane. Requiring 



-it) 



AttG 
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maximization of the proper entropy leads to the relation y/gooT = const. Thus, the proper 
temperature measured by a local observer will grow with the warp factor as one moves away 
from the visible brane. This growth prevents the flow of heat between regions of differing 
gravitational potential. We will further assume that the temperature is low enough so that 
the modes of interest are not strongly interacting. 



Let us first consider the field of interest to be, for now, a scalar, whose Hamiltonian will 



be 



2 ^ 



We then perform the usual Kaluza-Klein decomposition 



</.(x,y)=E0n(x)^, (35) 



where ipniv) satisfy 



and the inner product is defined as 

J'^ %V'n(l/)V'm(l/)e-'^«(^) = Sn^^. (37) 

The eigenmodes are linear combinations of Bessel functions, and were computed in p3[. The 



relevant components of the stress energy tensor are 

Tu = E ^My)'TS + ^0^(x)e^-^°(^)j-( ^"^^?5^^^ e''''"^^^^: 



= E l0^(x)e-o(4(M4^e--o(.)) + l0^(x)MM!. 
V 2 dy bo 2 bo 

Tqq and are the four dimensional energy density and pressure for the nth mode. In these 
expressions we have dropped mixing terms which will not contribute to the thermal average 
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and have used the free field equations of motion to simphfy T55. T05, which would give rise 
to a net fiow onto the branes, has vanishing thermal expectation value. Consider now the 
contribution from one, say the lowest, Kaluza-Klein mode n = in the limit of vanishing 
bulk mass (a fine tuned case when the bulk is not supersymmetric) . In this limit, the lowest 
mode is massless with a constant wave function. We find 

PBiy,t)=e'^^^y^i^',{y)^T\ (38) 

where again T is the proper temperature, as measured by a visible brane observer. Note 
that, as one would expect, the local energy density is weighted by the wave function and now 
has a profile which will only be consistent with our solution if it is trivial. The normalized 
wave function is given by 'ipQiy) = kbo{l — e~'^^''°), and we see from eq.(^) that, up to 
exponentially suppressed terms, the zero mode contributes to the energy density exactly 
like a field trapped to the brane, with its normalization set by the AdS curvature scale. 
Thus, the solutions discussed in the previous sections are valid for a massless bulk field in 
the limit where the temperature is too small to excite the higher Kaluza-Klein modes. It is 
straightforward in this case to see that the bulk pressure in the direction transverse to the 
branes vanishes and that bulk pressure just reduces to the four dimensional pressure in the 
same way that the energy density does. 

Before going on to the more interesting massive case, we point out that a massless 
bulk field is interesting in that it leads to no deviations from standard cosmology, which 
seems rather surprising at first. The lowest mode of a massless bulk field generates no bulk 
dynamics at all. Furthermore, its contribution to the expansion rate is exactly the same as 
in the standard cosmology. Contrast this to a massless field trapped to the brane which 
has non-linear contributions to the expansion rate, and which generates a shift in the lapse 
function 0. 
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V. THE MASSIVE BULK FIELDS 



Thus, we would like to extend our previous results to the case of a massive mode. 
Again, to make the problem tractable, we will assume that we are only exciting the lowest 
eigenmode. In principle it is possible to treat the sum of several modes without conceptual 
change. We will find solutions to Einstein's equations by perturbing around the massless 
solution via an expansion in the small parameter rriB/k. We will further assume that the 
bulk energy density is small compared to the brane tension, as we did in the previous 
section. This will allow us to drop time derivatives since they are higher order in p. Writing 
ipn{y) = C(l +"^B<^(z/)) leads to 

2 

Ttl 

6"{y) - 4kbo6'{y) = 6^(1 - -^e'^'^y), (39) 
whose normalized solutions are given by 

^(y) = kbo (^1 + ^(3 + 2e2^''o^ - e2'=f'o(22/-i)(i _ ^~2kbo^ _ 4^^^(^ ^ ^^^^ j ^^q^ 

where mo is the massive of the lowest lying mode, is the bulk mass and exponentially 
suppressed terms have been dropped. We impose the boundary conditions 'ip'^O) = V^'(l) = 0, 
since we have not included any interactions of this field on the brane. This assumption is 
rather unnatural since quantum effects will always induce operators localized to the brane 
pl| , |22|| . However, it is not difficult to modify this analysis to include these effects, which only 
complicate the algebra and do not lead to any conceptual changes. Imposing the boundary 
conditions we find that the mass of the mode is given by rriQ = — e^^^^"). This case is 
again a fine tuned case since the natural scale for the bulk mass is order Mpi, but is sufficient 
for our purposes. The average energy density is again just the canonical four dimensional 
energy density for a thermal fluid, as a consequence of the topological constraint arising 
from (|3|). 

Since the bulk profile is now no longer trivial, we will have to generalize our naive ansatz. 
For later convenience we choose the slightly different form 
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Where 

aiy^t) = Ao{y) + 5a{y,t) 
P{y,t)=Ao{y) + 6f3{y,t) 
b{y,t) = bo + 6b{y,t). 

All the variations are of order 0{^y^), so that their time derivatives will be suppressed. 
We may neglect the brane matter since their effects are controlled by a different expansion 
parameter, and the linearized results from the previous section may be carried over directly. 

It has been pointed out in that the Einstein equations simplify quite a bit, if we 
choose to work in terms of the variables 

V = Sp'iy, t) - {A',{y)/-^^ aiy, t) = Sa'{y, t) - Sp'{y, t), (42) 

which are invariant under gauge transformations which leave the branes invariant. In terms 
of these variables, the Gqo ^^nd G^^ Einstein equations are 

4AMr]{y,t)-r]'{y,t) = —K^5Too{y,t), (43) 

A',{y){Ar]iy,t) + a{y,t)) = ^'6n,{y,t), (44) 

and the combination Gqo + ble'^^^^y^Ga is given by 

^'{y,t) - ^A',{y)a{y,t) = -kV^»(^)6^(5Too(i/, t) + ^^^). (45) 

ao[t) 

We note that the average of the first two equations vanishes as a consequence of the topo- 
logical constraint imposed by the fact that we are working on a compact manifold. This is a 
manifestation of the point that to linear order in the perturbations there are no corrections 
to the averaged Einstein equations [0. We have the further simplification that both the 
energy density and bulk pressure have the same y dependence so that the solutions to the 
first two equations are identical up to an overall factor. We find 
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V{y,t) = -^^P4(t) (- {e^'^'^'^-'y^kboy (e'=^o(i+^) + 2cosh(A;6o (y - 1))) (coth(6oA;) - 1)) + 
g6ofc+46ofcy (coth(6oA;) - 1) smh(6oA; {y - 1)) sinh(6o%) + e^^''''^ (2 + e^''^'^) bokcsch{bokf smh(6o%)) , 

and 

a{y,t) = -3{l + ^My,t). (46) 

In deriving these expressions we have dropped the contribution from the derivative turns 
in the expression for the stress energy, as they only contribute at higher orders in the mass 
expansion. We may now check to see if the G55 equation is satisfied. This is the equation 
which accounts for the induced radion potential and leads to, in the case of brane matter, 
the unphysical constraint p = —e~^^^°Pi,. To the order we are working T55 vanishes and the 
constraint equation leads to the condition 

P4(T)=3P4(T), (47) 

which is of course not possible for a massive thermal fluid. Thus, we see that, at least with no 
external stabilization mechanism, there are no physical solutions with a thermal bulk fluid. 
It is simple to include the effects of matter on the branes since their effect is controlled by 
a separate expansion parameter. Since the Einstein equations are local and the constraint 
on the brane matter is independent of the fifth dimension, there is no way to even fine tune 
the relation between the brane and bulk matter to get a solution. 



VI. INCLUSION OF A GW MECHANISM 

To determine the effects of a bulk thermal fluid, we now include an additional scalar ^, 
which is responsible for stabilizing the space. This stabilization is accomplished by including 
a bulk as well as brane potentials for ^. An exact solution to this system without matter 
was found in ||2^. Here we will perturb around this solution, and following p], choose the 
simple bulk potential V{^) = |m|^^ and brane potentials Vi{C,) = mi{^ — ViY- The solution 
to the unperturbed coupled Einstein scalar field equations are approximately given by 
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Uy) = voe-'"''' Aoiy) = kh,y - ^vl{l - e'^^""^), (48) 

where e ~ A large hierarchy, without fine-tuning, is obtained by choosing e"'^'"' = 

{vi/vqY \ The Goo and Ga Einstein equations are unchanged, except for the fact that 
gauge invariant variable is now 

r^{y,t) = 5a'iy,t) - A',{y)^-^^ - ^^'^{y)Say,t), (49) 

Oq 6 

where 6C, is the deviation induced by the mass of the bulk matter. Note that there are still 
no jumps on the branes. That is, in these Einstein equations there are no implied delta 
functions induced by the shift in the fields, as they have all canceled once the leading order 
equations have been used. 

The G55 equation is now given by 

AMiMy,t) + <y.t)) + y [ioWiM - ayn\y.t) + {^y)/-^^^ = o + o{m'jk'] 

(50) 

Given that we have already solved the other Einstein equations we can solve the entire 
system by making the simplifying assumption that the brane potentials are stiff, so that the 
fluctuations of ^ on the branes can be taken to be zero. We may then use residual gauge 
invariance to eliminate 6C, in the bulk as well. Doing so allows us to calculate the shift in 
the radion due to the bulk fleld. We flnd 

Jo y bo ~ 24k'k^viboe'^^'^ ' '^^'^^''^ ISk^m^^ fcfeoM^^e-^^^o " ^^^> 
rrir is the mass of the radion which in the limit we are working is approximately given by 
mf. = ^k"^ {evoky e~'^^''^°~^^^ Note that the bulk fluid again behaves like a brane fluid, 

which couples conformally to the radion. 

The more interesting quantity is the change in the lapse function as this quantity tells 
us how the weak scale changes with time. Writing 

Mw{t)/Mw{to) ^ e-'/o '^"'(^■*)'^^ (52) 
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We find 

M^(t)/M;^(to) ~ 1 - ^2M2^I^y + ^;^) + (53) 

Notice that in the hmit where mo < T" < mi (which can be attained by fine tuning the 
bulk mass), and where our results are still valid, the contribution from the radion shift is 
subleading. 



VII. CONCLUSIONS AND OUTLOOK 

In this paper we have calculated perturbations away from the vacuum AdS^ compactified 
geometry due to the presence of a thermal bulk in the limit where T < tukk- We have 
shown that a bulk fiuid contributes to the expansion rate in a similar fashion to a fiuid 
which is confined to the visible brane and that the induced dynamics of the fifth dimension 
are controlled by the mass of the bulk field. The existence of the bulk fiuid does nothing 
to ameliorate the usual fine tuning between the brane matter when there is no stabilizing 
potential. In fact, without a stabilizing potential there are no physical solutions, at least 
that can be found by perturbing around the vacuum solution, when a thermal bulk fiuid is 
introduced. Upon introducing a stabilizing potential we found that the shift in the weak scale 
due to the existence of the bulk fiuid is parametrically suppressed relative to the contribution 
from a brane by m%/k'^. This stems from the fact that all the bulk dynamics induced by 
the bulk fiuid are controlled by this parameter. Indeed, we may conclude that the onset of 
dynamics in the fifth dimension, due to a bulk field, is triggered at a temperature near the 
mass of the first mode which has a non-trivial vacuum profile. 

Our results are only valid in the limit where T < rrij^j^, where ruj^j^ is the mass of the 
second Kaluza Klein excitation, so that contributions from higher excitations are suppressed. 
It would be interesting to find solutions where the mass of the first mode approaches the 
TeV scale, so that there is no need for fine tuning. This would entail solving the Einstein 
equations with the full Bessel function wave function source. Furthermore, one could then 
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also include more Kaluza-Klein excitations. This would allow for a study at temperatures 
high enough that the free energy begins to scale like T^. Presumably, as the temperature 
increases further, the local energy density near the Planck brane will approach the Planck 
scale (this can't be seen from our results since our calculation of the energy density was only 
valid for smaller temperatures where we could do perturbation theory). Once such energy 
densities are reached, black-hole formation seems inevitable and the change in geometry to 
AdS-Schwarzschild marks a phase transitions in the dual conformal field theory. 

ACKNOWLEDGMENTS 

The author is indebted to Walter Goldberger, Ben Grinstein, Rich Holman and Ted 
Jacobson for helpful discussions. This work was supported in part by the Department of 
Energy under grant numbers DOE-ER-40682-143. 



18 



REFERENCES 



[1] 
[2] 



[3] 
[4] 



[6 



[7; 



[9 

[10, 

[11 

\12 



L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) ||hep-ph/990522T| . 

P. Binetruy, C. Deffayet and D. Langlois, Nucl. Phys. B 565, 269 (2000) |Kepj 

th/9905012|| . 



N. Kaloper, Phys. Rev. D 60, 123506 (1999) ||hep-th/9905T[r: 



T. Nihei, Phys. Lett. B 465, 81 (1999) [|hep-ph/9905487| 
C. Csaki, M. 
ph/9906513|| . 



J. M. Cline, 
ph/9906523|| . 



Graesser, C. Kolda and J. Terning, Phys. Lett. B 462, 34 (1999) 
C. Grojean and G. Servant, Phys. Rev. Lett. 83, 4245 (1999) ^egj 



C. Csaki, M. Graesser, L. Randall and J. Terning, Phys. Rev. D 62, 045015 (2000) 



||hep-ph/9911i06 



J. M. Cline and H. Firouzjahi, Phys. Lett. B 495, 271 (2000) ||hep-th/ 000818^ 



R. N. Mohapatra, A. Perez-Lorenzana and C. A. de Sousa Pires, Int. J. Mod. Phys. A 
16, 1431 (2001) ||hep-ph/ 0003328 1. 



B. Kyae, |hep-th/00lll45| . 

J. E. Kim and B. Kyae, Phys. Lett. B 486, 165 (2000) ||hep-th/0005T39 |. 
J. Lesgourgues, S. Pastor, M. Peloso and L. Sorbo, Phys. Lett. B 489, 411 (2000) 
||hep-ph/ 0001086 



[13] B. Grinstein, D. R. Nolte and W. Skiba, |hep-th/ 00 12202 . 
[14] P. Kanti, K. A. Olive and M. Pospelov, Phys. Rev. D 62, 126004 (2000) |E§pj 
I ph/0005T46| . 

[15] M. Cvetic, H. Lu, C.N. Pope, |tiep-th/ 0002054) ; R. Ahendorfer, J. Bagger, D. 
Nemeschansky, |hep-th/0003117| ; T. Gherghetta, A. Pomarol, [liep-ph/0003129| ; N. 
Alonso-Alberca, P. Meessen, T. Ortin, Phys. Lett. B482 (2000) 400, piep-th/000"3248 



A. Falkowski, Z. Lalak, S. Pokorski, [hep-th/0004093| ; E. Bergshoeff, R. Kallosh, A. 
Van Proeyen, JHEP 0010 (2000) 033, |hep-th/0007044| ; M. Zucker, |hep-th/000'9083 



H. Nishino, S. Rajpoot, [hep-th/0011066| . M.A. Luty, R. Sundrum, [hep-th/0012158 



J. Bagger, D. Nemeschansky, R.-J. Zhang, tiep-th/0012163| . M.A. Luty and R. Sun- 
drum ^ep^thToTolTS^ 

[16] For a complete list of references see H. Davoudiasl, J. L. Hewett and T. G. Rizzo, Phys. 
Rev. D 63 (2001) 075004 ||hep-ph/0006im|] . 

[17] S. S. Gubser, Phys. Rev. D 63, 084017 (2001) [|hep-th/9912001 . 



R. Rattazzi and A. Zaffaroni, JHEP 0104, 021 (2001) ||hep-th/001'2'245 



N. Arkani-Hamed, M. Porrati and L. Randall, |hep-th/0012T48 . 

M. Perez- Victoria, |hep-th/010504^ . 
[18] I. Brevik, K. A. Milton, S. Nojiri and S. D. Odintsov, |hep-th/0010205 . 
[19] W. D. Goldberger and M. B. Wise, Phys. Rev. Lett. 83 (1999) 4922 [|hep-ph/9907447 



[20] P. Kanti, K. A. Olive and M. Pospelov, Phys. Lett. B 481, 386 (2000) [|hep-ph/0002229 



[21] H. Georgi, A. K. Grant and G. Hailu, Phys. Lett. B 506, 207 (2001) [|hep-ph/0012379 
[22] W. D. Goldberger and M. B. Wise, |hep-th/0104r7q 



[23] W. D. Goldberger and M. B. Wise, Phys. Rev. D 60, 107505 (1999) ||hep-ph/99072T8 



[24] O. De Wolfe, D. Z. Freedman, S. S. Gubser and A. Karch, Phys. Rev. D 62, 046008 

(2000) ||hep-th/ 9909 131 . 
[25] W. D. Goldberger and M. B. Wise, Phys. Lett. B 475, 275 (2000) ||hep-ph/9911457|] . 



19 



C. Charmousis, R. Gregory and V. A. Rubakov, Phys. Rev. D 62, 067505 (2000) jE^ 



th/991216q| . C. Csaki, M. L. Graesser and G. D. Kribs, Phys. Rev. D 63, 065002 (2001) 



| hep-th/0008T5T 



20 



